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in the development of computers used in fire control systems, 
there arises the requirement that two vector quantities (represented 
as scaler signals) be added vectorally at right angles. A method that 
was suggested to the author by Dr. Y. J. Liu was to use a closed loop 
system involving the squaring property of electro-magnetic microsyn 
units operating on a shaft. Although the instrumentation of such a 
system is fairly straight forward, the dynamic analysis of the system 
is made complex because of the non-linearities involved. This thesis 
will be in the main a study of the dynamics of the proposed system. 

After a brief discussion of the theory and operation of the 
microsyn units to be incorporated in this system, the performance 
equations along with stability conditions for an ideal system are 
developed in Chapter II. Chapter III then develops a generalized 
approach to the design of a practical system. То keep the dynamic 
qualities of the system constant a square root generator has to be 


incorporated in the feedback loop. What was done in Chapter III was 
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to develop a method specifying how close to an ideal square root 
characteristic we have to be in order to stay within pre-assigned 
tolerances of the natural undamped frequency and the damping ratio 
associated with the system. Chapters IV and V present transient and 
steady state analyses similar to the types used in linear second order 
systems. The final chapter shows the complete design and results of 
a laboratory test system using the general design criteria of the 


previous chapters. 


Thesis supervisor: Dr. Walter Wrigley 


Associate Professor of Aeronautical Engineering 
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CHAPTER I 


Introduction to the Problem 


In the instrumentation of systems related to a wide variety 
of problems, there often arises the need to add vectorally two quan- 
tities at right angles. In general, the two vector quantities are 
represented as scalar signals, and it is the requirement of some 
instrument or system to treat these signals as vectors and perform 
the right angle addition. In this thesis, we shall represent the 


two signals as voltage quantities E, and Еу. The desired result will 


then be Eo = {ερ + E^ where again Eg is 8 scalar quantity. It is 


to be noted here that we are not interested in the actual position with 
respect to inertial space of the two vectors, but merely in restoring 
their right angle relationship with respect to each other. 

The quantities E, and E, will fluctuate either periodically 
or haphazardly, but they will never change direction. The reason why 
E, and Еу can only be fluctuating direct voltages of the same sign will 
be shown later in the chapter. A major portion of this thesis will be 
devoted to a study and analysis of the dynamics associated with the 
proposed system. In short, we shall attempt to answer the two questions, 
"how does the system respond to fluctuating inputs?"- and given certain 
tolerances on how well the system should perform, "how can we design a 
practical system?" 

One method in use at present to obtain the result 


nos (52 + Еу“ is to use a resolver nulling servo combination. The 


system is shown schematically in Fig. 1-2. 





DESIRED RESULT. 


Е ж 
j E 2 2 


Е, 
FIGURE 171 
RIGHT ANGLE RELATIONSHIPS OF QUANTITIES TO BE 
COMPUTED 







SERVO 
MOTOR 


RESOLVER SHAFT 
q == E = wif) ) 






AMPLIFIER 


FIGURE 1-2 


USUAL METHOD OF SOLVING PROBLEM 


NULLING SERVO 





Е, апа Ey are root mean square values of alternating sinus- 


oidal voltages appearing across the stator windings I and II. Eo] nd 


Ер аге voltages appearing across the rotor windings as shown. The 
axes of the two stator windings are located 90 mechanical degrees apart 
as are the axes of the two rotor windings. Because of the sinusoidal 
relationship of flux density in the air gap as a function of angular 
position of the stator periphery, the following relations hold for 


open circuit rotor windings. 


II 


Eol - E sin @ + E, cos O! 1-1 


Ес2 Ex cos € + Ey sin Q l-2 


In the case of the resolver nulling servo combination as 
shown in Fig. 1-8, the shaft of the resolver is driven in a direction 
by the servo motor so аз to null E,j- If the nulling г а amplifier 
has a large enough gain, Eo] Can, to a first approximation, be con- 
sidered zero but it is understood that Ej] = 9 only if the gain of 
the nulling servo amplifier is infinite. Thus equations 1-1 апд 1-2 
become 

Ео2 = Ex cos Q + Ey sin Q 125 
0 = - E sin O + Ey cos Q 1-3 


Squaring equation 1-2 and 1-3 and adding, we obtain 
E ЭС = 2 2 
0 1 s Еу 


and thus we have our desired result 


Boo 
2 2 
Eco z Bo = E + B; 1-4 
This system as snown in Fig. 1-2 1з а true servomechanism 
and its characteristics will be used later as a comparison showing the 


differences between this true servomechanism and the closed loop system 


to be described in this thesis. 
We shall hereby name the proposed system the orthogonal vec- 
tor component summing system. A simplified diagram is shown in Fig. 1-3. 
At this point, we shall neglect the damping and inertia terms 


and merely show how the steady state result is obtained. In Fig. 1-3 


the two inputs E, and ES &re acted upon by an operating component (oc) 
such that an output torque is produced on a shaft (torque summing 
member) which is proportional to the square of the input signals, The 
shaft rotates through an angle A rot) to which is attached an operating 
component such that an output voltage is produced proportional to the 
angle of rotation. The square root of this voltage is taken and fed 
to the shaft through a third torque producing device which depends 
on the square of its actuating input quantity. 
In the steady state, that is, assuming all transients have 
died out, we have the three torques М,, Mo, and М; acting on the shaft. 
M + M, я M, 1-5 
ог га 
2(0c)1 le25u] 62 + 5(ос)2 |в2м) в” = S(oc)3 le?;u| Eo ia 
Assuming the sensitivities giving torque out for voltage 


square in are equal, we arrive at the desired result, namely, 
2 2 
Eol х Еу + Ey 


Е = | Жы + E, 1-7 


The square root function in the feedback loop cascades with 


the squaring device to produce & linear device in the feedback 


x 
Refer to the Appendix for explanation of terminology used. 
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loop. Equation 1-6 shows the very important result that in the steady 
state, the correct answer to the problem will be given regardless of 
the neture of the device used in the feedback loop. However, in the 
next chapter where the complete performance equations of the system 
will be developed, we shall show that a perfect square root function 

is required in the feedback loop if the dynamic qualities of the system 
are to remain constant. The remainder of this chapter will be devoted 
to a description and theory of the torque producing generator and the 
angle input voltage output component commonly called the signal 
generator. 

Both the torque generator an? signal generator belong to a 
family of devices known as microsyn electromagnetic s These 
units consist of a laminated ferromagnetic rotor of 180° symmetry with 
two poles and a laminated ferromagnetic stator of 90° symnetry having 
four poles about which are wound various arrangements of electrical 
coils depending upon the application of the unit. For the types of 
microsyns used in connection with this thesis, the rotor is allowed to 
rotate * 109 from & zero rotor position, which is geometrically defined 
when a rotor pole is half way between two adjacent stator poles. Fig. 
2-4, shows the general construction of a microsyn. torcue generator. 

The torque generator produces a torque which, in the ideal 
case, is proportional to the square of the current through one axis 
and is independent of rotor angular position over the useful range. 

Meg) = S(eg)[s2sm] И 54 

We shall attempt to clarify this relationship from energy 

considerations in a magnetic field. 


As a starting point, we note the fact that the total 
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electrical energy into the system equals the magnetically stored 


energy plus the work done by the system, neglecting heating losses. 


That is, 
(En) (oec) 3 (Еһ) (пав) μα ш 
The instantaneous electric power delivered to the torcue 
producing coil is ei. Thus, after a time t, the electric energy 
into the system, is 
t 
(En) (e1ec) = / ei dt 
According to Faraday's law 
ena) = - NSP 1-10 
where 
elina) 3 voltage induced by motion 
N = the number of conductors 
ф = magnetic flux cut by the conductors 
Hence be 
CERIS -S'n ag idt -/ № аф 1-11 
% $ 
Where 
0, iscfüux 3t b 
P> is flux at t = ty 


We shall go from Ø} to $2 by two paths as shown in 
Fig. 14. First, following path 1, the current is maintained constant 
and hence 
(En) (elec) ^ Ni (95-99 1-1; 
This will correspond to a certain rotation of the rotor 
to a point corresponding to the flux Φ 2: Next, following path 2, 


we let the current drop to zero, and if we neglect any hysteréses 





FIGURE l-5 


AMPERES VS. FLUX IN A MAGNETIC MEDIUM 





effects, the flux too will drop to zero. At this point the rotor is 
rotated through the same angle as in the first case. However, no work 
is done because there is no force. Then the current is increased until 
the point (Pp, ii) is reached. In the second case, the magnetic energy 
can be written in the following form. 


-/ маф + f лав 


Ф 


+ JL - y^ 1-13 


о 


(En) (elec) 


1 : 
Going from ф 1 to zero, denote Kj - ф where Kj is the 
1 
slope of this line. Similarly, going from O to $5, denote K5 Ξ D 


the slope of this straight path. Then 1-13 becomes 


2 ф 
(En) (elec) = ку / фаф- NK IS 1-14 
о 


ο 
But 
Therefore | 
(En) - wil аф -νὮ e 
(elec) Pl М и 
N il N il 2 
Е 2% ka 
EXIT 
(En) (e1ec) то (Φι - ф2) 
or finally 


(Ел) (елес) - : i ($5-$) те 


10 





x 


Equation 1-12 represents the area under path 1 while 1-15 
represents the difference between the areas under the two sections of 
path 2. 

Since 1-15 is exactly half that found in 1-12, we conclude 
thet the mechanical work done and the magnetic energy change are equal 
in magnitude for а given amount of input electric energy. 

The fact that the magnetic energy stored in the magnetic 
field equals the work done under the conditions previously described 
enables us to calculate the work done by finding its equivalent in 
terms of the magnetic energy storage. Then, the torque can be found 
by taking the derivative of this amount of energy with respect to the 
rotor angle. 

Before proceeding further, we shall make the following 
assumptions: 
l. The reluctance of the magnetic circuit is concentrated almost 
entirely (over 90$) in the air gaps. As a consequence of this assumption, 
we shall assume that the entire magnetic drop around the magnetic loop 
produced by the sources of magnetomotive forces at opposite poles occur 
essentially across the air gap. Thus, we can assume to a good approxi- 
mation that all the magnetic energy is store? in the air gap. 
2. Leakage flux will be neglected. 

From reference (2), the magnetic energy stored in a given 


volume is given by 


EDI 2 
(En) (mag) za LIT Г Нав 1-16 


1 
where V is the volume, B is the flux density. As mentioned above, 


this energy will be assumed to be stored in the air gap. 


Hence: 


11 





В 
ү ГА 2 ү 2 2 
(En) (mag) = 17 Ue Вав = атр (Во - By) | 


В 
1 
where із the permitivity of air in MKS units. 
2 2 
LA A. : о "ντ ο - $4) 
E = | uoo к= 2 1 
(En) (mag) ЛИ A ενα; 


1-18 
where L is the distance between the rotor and stator. 
A is total area between the rotor and the two opposite 


stators. Since 


б - fo c. pe 
τ reluctance E (mag) 
L 1 


2 
Е = —F T (change in mmf) 
(En) (mag) ӨЛ ДА 8612) 


2 
oA U(mag)^ — 1-19 


8 TL 


where U(mag) is the magnetic drop across the air gaps. This drop 
equals to a good approximation the impressed magnetomotive forces. 


Since mmf = ДЖ М№і, ме have finally 


22 
= 2T uL AN i ` = T= - 
(Επ) (πας) = Т W Neno. turns 1-20 
The torque M then is Ow where A(rot) is the rotor angle. 
A( rot) 


Thus it is shown that the torque is proportional to the 
current squared. An important consequence of 1-20 is the torque is 
independent of the direction of the current. Thus alternating current 
can be used as well as direct current to provide torques. This also 
shows why only absolute values of Ey and Ey can be considered in the 


orthogonal vector summing system. 


12 


The change of area A has a linear relationship to a change 


in rotor angle A( pot): Thus 1-20 becomes 


С 2п Мо КАМ i - -27 Ho N^i^K A(rot)_ 
E Г 


(En) (mag) L 


where K can be calculated from the dimensions of the particular 


microsyn. 
Finally, 
оса 2 
" =- 2THgKN i = с ‚2 1 1-21 
Κερ) = ОКБ L (ta) [3°30 ] 


where we shall define 


2Ί µο ΚΝ΄ 
ООЛ 0 


The signal generator provides an open circuit voltage which 
is proportional to the amount of angular rotation of the rotor from 


its null position. Thus 
(зе) = 5(зв)| А;е | А (тоф) 1-22 


Fig. 1-6 shows the details of the signal generator. There 
are two windings per stator pole anc these windings are connected in 
two series circuits, each circuit using a similar set of windings. 
The primary circuit is excited by an alternating current which sets 
up flux patterns as shown. If the previously mentioned torque generator 
has equal excitation currents on both axes, the flux pattern will be 
the same as that set up by the primary excitation current of the sig- 
nal generator anc thus the resultant torque on the signal generator 
is zero since the torques on the two axes are equal and opposite. There 
will be no torque contribution from the seconcary winding since this 


winding will carry no current. The series windings of the secondary 


13 
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MICROSYN SIGNAL GENERATOR 
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coils are so arranged that the instantaneous induced voltages in these 
secondaries are added in opposite phase for adjacent poles. 
Thus, 

uc NE NEL Uo О 1-23 
where Eg 3s the rms induced voltage in coll n. n2 1,2, 3, or 4 

When the rotor is in its null position (i.e. midway between 
two adjacent stator poles), we see from 1-23 that the signal generator 
output voltage is zero because all the in?uced secondary voltages are 
zero due to the symmetry of the system. 

A positive rotation of the rotor will be defined as a 
rotation tending to increase the overlap area between the rotor and 
stator poles 2 and 4 and correspondingly, to decrease the overall over- 
lap between the rotor and poles 1 and 3, as shown on Fig. 1-6, Thus, 
fcr a positive rotation, the reluctance of the magnetic circuit inclu- 
ding poles 1 and 3 increases while that of 2 and 4 decreases. Hence, 
since the magnetomotive force at all the poles is constant, the flux 
in poles l and 3 decreases while that of 2 and L increases. It is to 
be noted that the change in reluctance and therefore the change in flux 
is a linear function of rotor angle. The incuced root mean square 
voltages are proportional to the flux for constant excitation frequency. 
Thus for a positive rotation, the induced voltages in poles 2 and 4 
ere the larger and the output voltage has a phase similar to the phase 
of the voltages in poles 2 and 4. For negative rotations, e(oyt) 18 
of opposite phase and of magnitude given by 1-23. 

The flux in the poles of the signal generator have very 
nearly the same phase as that of the alternating excitetion current, 


the only difference in phase being that attributable to the core losses. 
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This flux is directly proportional to the number of primary coil turns, 
and the magnitude of the excitation current and inversely proportional 
to the circuit reluctance. The secondary induced voltsge has a 90° 
phase shift with respect to the flux and its value is proportional to 
the magnitude of the flux, number of secondary turns, and excitation 


frequency. We can thus define the output voltage as 


e(out) = S(sg) [A5i,n;e] l(ex) n(ex) А (го%) 1-24 


2(sg) [А ;e | A rot) 
where 


S(sg) |4;1,n;e] 1(ех) (ех) 7 °(se) [Ase] 


when the excitation current and frequency are fixed, and A(yot) is the 


angle of the rotor measured from the null position with a positive 


rotation having been already defined as a rotation which increases the 


area under poles 2 and 4. 
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CHAPTER II 


Proposed Ideal System 


The orthogonal vector summing system to be described in this 
chapter will be considered ideal for the following reasons: First, the 
square root generator will be considered as producing an exact square 
root over the entire operating range. Second, all uncertainty torques 
such as those attributable to friction, hysteresis, leakage flux, plus 
interfering torques such as the acceleration of the case with respect 
to inertial space will be neglected. 

Fig. 2-1 shows a line schematic diagram of the proposed 
assembly and Fig. 2-2 shows a functional diagram of the orthogonal 
vector summing system. Since this thesis is intended to be primarily 
a study of the dynamics associated with the proposed system, the des- 
cription of the actual assembly will be mainly to establish a clear 
picture of what physical processes are actually taking place. The 
analysis to follow will be largely mathematical, but nevertheless the 
physical interpretation of the derived performance equations will be 
continually stressed. 

The assembly shown in line schematic form in Fig. 2-1 is 
mounted inside a sealed cyclindrical container and designed to operate 
іп an upright position. The torque summing shaft contains the rotors 
of the microsyn units as well as the rotor of the viscous damper. The 
fluid in addition to providing a damping torque, also serves to reduce 
friction torques at the pivots by buoying up the shaft by a sufficient 


amount to remove the weight of the shaft from the pivots. A float is 


17 


№315 45 
 ΝΙΝΝΗΣ Ч80193Л IYNOJOHIYO JHL 30 WVYOVIG ОПУКЗН)5 JNN 


|-ό 38 Πν η 


YOLVY ANID 
1УМ915 








чаамуа 
SNOOSIA 


SYOLVY ANID 
SINOYOL 


18 


Н11545 ӘМІКНП5 YOLIZA WNOXOHIYO 1V30| 30 NVUSVIO TWNOILONAS 


93 





295 
[9:υ]{65}ς 
Qon 
(65) 
| 
 ——— а 


indui «ρου δωος 10500 


0-0 180913 


V I -— 
| (vo)[(09)-1],, 
(Уу Y 


VILYIANI 





Є 63 


07 


(WS } ) 














7 гә) 
“(324006 1М324400) 
Aq п*әуїзә) 
*(394n0S LNIYYNI) 
ε | 
1.9] (59) 
τ [1:9] S 





(324n00S ΙΝ388Π9) 


?3 


19 


clamped to the shaft to provide the necessary amourt of lift. The 
damper consists of two concentric cylindrical surfaces one of which 
rotates with the shaft while the other surface is fixed to the case. 

А damper adjustment is provided which changes the amount of surface 
between the concentric cylinders. Because the viscosity of the damping 
fluid is critically dependent upon temperature, the unit is maintained 
at a constant temperature plus or minus a certain tolerance by means of 
electrically heated coils. A thermostat is provided which sets up a 
cycling process in the application of electric power to the heating 
coils. А portion of Chapter V will be devoted to an analysis of how 
changes in damping effects the dynamics of the system, and thus 
determine criteria for specifying how well the temperature must be 
regulsted. 

We shall now develop in detail the performance equations for 
the ideal system referring to Fig. 2-2. Because of the fact that we 
shall regard the input signals as voltages, a current source has been 
provided to drive each of the torque generators. 

By Newton's second law, the product of the moment of inertia 
of the shaft plus rotors times the angular acceleration of the shaft 


eguels the sum of the torques applied to the shaft. 


(А [-(са)] (04) * Arot)) 7 Ca(vis) A(rot) * M(tg), * M(tg)2 


* M(tg)4 2-1 


Where I = moment of inertia of shaft including rotors of the damper 
and microsyn units. 
A(rot) = angle through which a rotor has turned from its null 


position. We shall assume all rotors are lined up to 
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have the same null angle with respect to case coordinates. 
A positive rotation of A(rot) has been defined as that 
direction tending to increase the overlap area between 
the rotor and stator poles 2 and 4. 

А[т-(са)] (OA) = the component of angular acceleration of the case with 
respect to inertial space about the output axis, that 
is,about the direction indicated by a line parallel 
to the shaft. 

Са = The damping coefficient of the system consisting of 
the following two parts: Cg@(fixeq)» the damping co- 
efficient of the fixed section of the viscous damper 
and Cq(var) the variable section of the viscous 
damper. 

Torque generator 1 and torque generator 2 will be actuated 

by the input signals and will be defined as producing positive torques, 

i.e. torques tending to produce a positive rotation of shaft angle. 

On the other hand, torque generator 3 will be actuated by the output 

of the system and will be designed to produce a negative torque. 


Rearranging 2-1, 


І (гоь) * Ca A(rot) 7 "(tg)3 - M(tg)i * "(tg)2 "° ` ШЕСІ (о) ** 
From chapter I, we have defined 
M(tg) = S(tg)[i°3™] 15 $25 
Define 5(сѕ) [е;і | = i sensitivity of the current source 
for voltage input current output. Hence 
M(tg) = (tg) [22 5м} S^ (es) [e; E eas 
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_ Бо 


А 1 - πα. 
Define s 6) per , Sensitivity of square root 


1 
generator for (voltage)? in, voltage out. 


From Chapter I, we have defined the signal generator sensi- 


Es 
tivity as 5(зв) [А;е | = ES where 


°(sg)| Ase | й "(sg) [A5i,n,e] 1(ех) П(ех) 


where і (ех) is the excitation current of the signal generator primary 
winding and П(ех) is the frequency of the excitation current. 
Using the above definitions for the parameters of the system, 


2-2 becomes 


! 2 ра 
D ос" (дз [зн] ? (es)4|e;i] cine 
2 2 
S(tg)3 [5m] SC Ἐς б 3(tg)2 [1^;« P (с)2|еҙ1 | 5 
a5 А[т-(са)] (ОА) 2-5 
біпсе 
"(491 7 52 = (ы); 7 (tg) 
апа 
(cs) T Sea T a 
I A(rot) * Ca A(rot) + Sg) [155] S (es)[esi] EA 2 
Sea ыы Ы : 83 
- I ЖЕСІ (OA) 2-6 
А150, 
во = 5" (sr) [eà;e] S(sg) [A;e] ^(rot) = 


22 


Using 2-7 and defining 
- 2 2 aes 2-8 
S(ovss) | A;M| 7 S(tg) [i> ,M| (ез еҙ) - (sr)|e?;e| (sg)|A;el 
the sensitivity of the orthogonal vector summing system for shaft 


angle input torque output, 


2-6 becomes finally 


2 
I on) + Cy A (rot) + S (ovss) |A;M] (rot) Ё S(tz)|[i2;M]” (es) |e;i | 
(2.2). ГА кы ον 
2-9 


In the discussion of the ideal system to follow, we shall 
assume that the acceleration of the case with respect to inertial 
space is zero. This acceleration produces an uncertainty torque 


which will alter the desired result. Thus, equation 2-9 simplifies to 


i Á( rot) + Са ο) + S(ovss)|A;M] A( rot) = (ъв) [1574] Ses) [ess] 


СЕ Я 5,2) 2-10 


Relating A(pot) and ET and the derivatives, we obtain 


Eo” 


EC 


S" (sr) [е®;е] (sa) [45e] 


φ 
— 


πα 
— — ο ο αν ο κ шз 


S (ar) е ;е| (за) [А;е | 


А (тоб) 


ο. 


where the independence of the sensitivities of the signal generator 
and square root generator with time are assume for the range of fre- 


quencies to be considered in this thesis. This statement will be 


23 


clarified later. 


Using the above expressions for A( rot) and its derivatives, 


2 2 а ee 
T — Eo* + Са d. Eo + S(ovss) [asm] Eo = 


дї dt 
S(ovss) [A5] E + Ey?) 2-12 

ог 
І ο. 52. Са а pts 82-562. By? 2215 


м 


S(ovss) [A;M] dt S(ovss)[A;M] dt 


Equation 2-13 is a linear differential equation of the 
second order with constant coefficients in the ideal case. Using 


the standard form of the second order linear equation, we have 


Ñ 2. Е2 2 Е2 Σ 
1 45 ұ2, 2(08)(0уѕзз) d g^,E? = Е + Е 2-14 
о + О 
Wn? dt Wn dt 
where W. = the natural radian frequency of the system. 


CPN uss) the ratio of the amount of damping present in 


the system to that amount required for critical 
damping. 
(DR) and Wn, When known for a second order system, completely determine 


its dynamic behavior. Comparing 2-13 and 2-14, we obtain, 


I 


1 и 
Wy? S(ovss) [A ju] 
iu fu [З(оузз) [A38] 
п 1 2-15 


ог 
258) (ονςς) Ж ба 
an S(ovss) I 
or 


Ca | 
(DR) (ovss) = -F V T S(ovss) [uti] 2-16 
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Since the two input quantities Ey and E, react upon the 
shaft through identical operating components, we can let 
Е + К = Ein^ in further discussion of the dynamics of the system 


without any loss in generality. Thus, 2-13 becomes 
2 - 
1 q m 


КИН 22 25-222 5 — Ca ___ d. КО. + B = Es 2-17 
S(ovss)|A;M | dt п 


Э(оуз) [a;M | 


or its equivalent form 


2 2 DR) а 2 Pere 2 2-18 
Wop een) 2: << - B 
%,2 dtZ © Wn dt O © іп 


Where Ej, is understood to be a function of time. 


In the steady state, 2-10 becomes, 
2 2 2 
= IB E 2-1 
S(ovss)|A;M | A( rot) P (tg) [μήν] Š (cs) [е;4 (% ay, ? 
or 


S (er) [e;e] °(sa) [Ase | A rot.) - Еу : Es 2-20 


From the steady state value of 2-13 we obtain 
2 = y 2 2 S2 1 5 А 2-21 
Eo^ 7 Ey € B^ - 9 (sr) [e*;e| (sg) [Азе | (rot) 


Thus, it can be seen that Eo” and hence Eo, is not 
critically dependent upon the value of A(yot), If the product, 
S Cer) [eise] sg) [Ase] is changed, A(rot) merely readjusts itself, 
within its constraints of πο. to a new steady state value corres- 
ponding to the required steady state.value of К. However, it is 
very important to note at this point that a change in the product, 
S*(sr) [ее] S(sg) [ase] > does not alter the steady state conditions, 
the dynamics of the system are definitely altered as can be seen from 
2-13, 2-15 and 2-16. The fact that a perfect square root function can- 
not be generated means that the sensitivity 5. will be continually 


varying over the operating range and the next chapter will go into the 
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details of the relationship between the dynamics of the system and 

the nonlinear device. It is to be noted here that the torque generator 
can be regarded as an almost perfect squaring device for the range of 
actuating currents used here. That is, the B-H curve is linear for 
the range of currents used, (0 to 40 ma) but there are small uncer- 
tainty torques due mainly to hysterises and to reaction torques. 
Reaction torques are present when the torque produced is not completely 
independent of rotor angle. For the present, uncertainty torques will 
be neglected. 

When the sensitivity product S^(sr) [ed ;e] "(sg) [A;e ] is 
made very high, A(.9,) decreases to e point where it becomes essentially 
an error signal analogous to that of a conventional servomechanism, 
similar to that shown in Chapter I. The system will be analyzed and 
designed, however, from the viewpoint that A(.95,)can be made to cover 
the full 10°. This condition is what makes this system unique from 
the conventional servo systems. 

A very important point on the stability of the system must 
now be brought up. At the null position, a condition of instability 
exists such that if A( rot) goes the least bit negative, the system 
drives itself into the stops provided just above the 10° mark of 
rotation. 

Fig. 2-4 shows š stable condition where the rotor has 
rotated through a positive angle. The rotor hes rotated to 4 position 
so that the signal which it procuces causes a current i, to be equi- 
valent to iin ano thus the positive torque produced by ijn equals the 
negative torque produced by ig. М А(роь) by some disturbance goes to 
A(rot) + 4 A(rot) with ijn the same, then ip goes to іо + d іо апа 


produces a torque tending to ру11 А(гоф) + Я A(rot) back to its initial 
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FIGURE 2-4 


TORQUE GENERATOR IN STABLE REGION 
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value of A(rot). If A(rot) is disturbed to a position A(rot) - d (του) 
ijp constant, then ig becomes ig - d igo and the torque produced in the 
1, 3 axis is reduced thus allowing A(rot) 7 а А (гоё) to rotate back to 
A(rot)- This again clearly shows the stable condition of Fig. 2-4. 

The configuration shown by Fig. 2-4 corresponds to a potential minimun 
which, of course, will thus yield a stable system. 

Fig. 2-5 shows the condition of instability. Assume, initially 
the conditions of Fig. 2-5. -Á( rot) produces a signal which will pro- 
Әсе the same current, io, а5 + A(rot): The condition as shown can 
exist just like a pendulum at the top of its arc. However, if 
- A(rot) is disturbed to a position, -(A(rot) + d Αρ, 

- A(rot) - 4 A(rot)» to goes to ip + d ig thus increasing the negative 
torque which in turn increases - A(pot) and the system will rotate 
negatively until finally hitting the stops. Оп the other hand, if 

- A(rot) is disturbed to - (A(rot) - 4 A(rot)) OF - A(rot) + 9 A(rot)> 
io goes to io - ἆ ip and the negative torque produced by i, is 
decreased thus tending to further reduce - A(yot) until the null is 
reached. A positive éngle is built up on the stable side until finally 
the conrition of Fig. 2-4 is reached. Here, we heve represented a 
potential maximum condition, it being similar to а physical pendulum 
at the top of its arc. The slightest disturbance in eitner cirection 
results in drastic changes in the motion of both the given computing 
system and the penculum. 

As & third case, assume initially there is no input excitation 
and DS = 0. Next assume e builds up linearly from zero. As iij 
builes up from 0 to d iin, there is a positive torcue pro-uced by the 
current flowing in coils 2 and 4 (see Fig. 2-4). At the same time, 


there is no current ip in coils 1 and 3 ane hence no negative torque 
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FIGURE 2-5 


THE UNSTABLE REGION OF SHAFT ROTATION 
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is produced. This assures a positive rotation and a stable condition. 


Finally, consider the case where A(yot) * O and ii, = O. 


If A(rot) is disturbed to a position - d A(yot), there results a signal 
producing an output current d ig and a negative torque which will in- 
crease the negative angle. Thus, we have the unstable case where A(rot) 
progressively increases until the stops are reached. 
We conclude from the above analysis that a constraint has 
to be applied to the shaft so that negative angles as shown in Fig. 2-4 
cannot be built up. If the constraint in the form of a mechanical 
stop is placed exactly at the null, then the system will produce an 
Ко for Ey and Ey even in the limit as Ek, and Еу both approach zero. 
If the stop, however, is placed at a point allowing a small negative 
angle - d A(rot), then the chances are that initially the rotor will 
have passed into this negative region in the absence of input signals 
as was pointed out in the last case of instability. Thus, it will 
take a positive torque produced by the input signals just barely 
greater than the negative torque produced by - а A( rot) in the feed- 
back loop to rotate the shaft into the stable region. It is difficult 
physically to put a stop exactly at the null, and in the last chapter 
it will be shown why a stop can be put on the positive sice of the 
shaft rotation and still be able to solve the orthogonal vector addition 
even as the input signals approach zero. This method will be feasible, 
however, only if the swing of A(yot) is fairly large in comparison to 
the angular position of the constraint, and also a certain amount of 
tolerance is allowed on the dynamical characteristics of the system. 
Equation 2-18, while linear when expressed in terms of the 
squeres of the input and output signals, is nonlinear when expressed 


in terms of the input and output signals themselves. 
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An insight into the complexity of the nonlinear form of 


2-18 can be gained if it is rewritten as follows. 


a 2 dEo 

Са - 25. 

dt ο ? at 
d^ к? = 2E, Бо.  2(9Ео)2 
dt? 7 412 dt 


Hence 2-18 becomes, upon rearranging, 


E, “о , (9.2. 2(08) м. Е № + Wn. EQ’ = "m Ejn (t) 222 
о 22 GEO 3 в. 


It is to be emphasized here that W, and (DR)refer only to 
the linear second order system defining the square of the actual input 
and output quantities. Since the linear equation defining Бо” 18 
convenient to use, practically all the information to be obtained 
from the system will be in terms of the squared functions. The last 
step will then be to analyze the square root of the response function 


associated with the linear second order equation. 
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CHAPTER III 


The Square Root Generator Problem 


The performance equations of Chapter II were developed 
assuming an ideal square root generator in the feedback loop. As we 
have already mentioned, it is very difficult to find in nature a phy- 
sical device that will generate such a curve as shown in Fig. 3-1, 
although of course, there do exist many non-linear devices or means 
of synthesizing such devices. 

Equation 2-8 defines the overall sensitivity of the system 
as 

B - 5 aloe T 54 ые” ; 
(ovss) [AM] (tg) [i ;M] (сз) [е;1| (sr) |e ; | (sg)| Ase| 

Since the torque generator follows a rather precise square 

curve for the range of currents used, only a perfect square root function 


will maintain S(ovss)|a;m | constant over the entire range. Physically, 


S(ovss)|A;M| is analogous to the spring constant K(sp) of a mechanical 
2 
harmonic oscillator. This is shown clearly in the equations for the 
natural frequency and damping ratio where 

S(ovss) fa;m] 

I 

Са I 
(28) (оузз) = 2 I S(ovss) [Ash] 


Wn 





where W, and (DRINS) refer to the performance equation relating the 
squares of the input and output functions. It is to be noted again here 
that the static output of the system is not dependent upon 3(оузз) [ АЗМ | 


If instead of a square root function we had 4 linear function, S(ovss) [^;X | 
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NUMERICAL CONSTANTS FOR 
THIS FUNCTION HAVE BEEN 
DERIVED IN CHAPTER YI. 
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FIGURE 3-1 


IDEAL NON-LINEAR CHARACTERISTICS IN FEEDBACK LOOP 
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would correspond to a very stiff spring for large rotor angles, and 
hence the W, of the system would be high, and the (DR) (ovss) would be 
low. For low rotor angles, the reverse conditions woulr be true. 

We thus have to resign ourselves to the fact that in an actual 
physical system Wp and (DR)(ovss) are not going to remain constant. Wy 
and (DR)(ovss) corresponding to the ideal system will now be defined 
as a reference natural frequency and a reference damping ratio, 
abbreviated, 

Wn(ref) 8nd Е, 

In attempting to solve the square root generator problem, 
we shall apply a method of perturbations to a simplified system and then 
critically examine the conclusions obtained from such a system. 

The simplified form of the perturbation theory will be a means 
of linearizing the non-linear equation relating Eg апа Ein. The pro- 
cedure will be as follows. We shall assume an arbitrary static operating 
point. The input signal will then be varied Бу а smail amount and the 
corresponding small change in the output determined. Expansions by 
Taylor's series will be used when necessary and becéuse the variations 
will be very small, A 's of higher огӣег than one will be neglected. 
When the so-called perturbed equation is subtractes from the original 
equation, a linear equation results. 

Fig. 3e2 shows the orthogonal vector summing system in some- 
what simplified form where instead of square root and signal generators 
we have merely specified a device whose output, E, is a function of the 
rotor angle, A. 

From Fig. 3-2, and using Newton's second law of motion 


ЕСКИ е zT 
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ог 


I Á + Ca hae Ыр =o My з=? 





о S(tg)| i2] Sees) е;1 | Eo" 
= τ; 52, үг .1 B, 2 
М = 5 (в) |2 | (сѕ) [е;і | іп 


We shall abbreviate the shaft angle A(rot) to A in this 
chapter because of the large number of times that it will be used. 

Using the fact that 

Е = f(A), 3-2 becomes 
ІА + Са А + S(tg) pem) [(4)] ИЕ S(tg) S4 Ces) E“ Cin) 355 

where we have assumed that the torque generator sensitivities and 
current source sensitivities ere equal. Using the method of pertur- 
betions, we let E;, assume a new value Ejn + A Ein where A Ein is a 
small change in the static value of the input signal. Then, corres- 
ponding to this change the output rotor changes to A + MA and the new 
output can be written f(A + A A). 


Equation 3-3 for this change in input then becomes 
2 А 2 
Ее ΔΑ)] = "^ (tg) 5 (св) Еш” АЕ: п | 


3-4 


-- + 9 + 9 52 
К ШЕЛЕК кид 008) 498) 


Expanding f(A + A A) in terms of f(A) by means of Taylor's 
series we obtain, 


f(A) 


d 
КА ДА) = f(A) + ЗА Αα Αν ο. 3-5 


We have neglected A's of higher order than one. 


Thus, 
[4+4]? = [#4] + 2 ГИДА... 3-6 
where -— _ ΤΡ 
аА 
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апа 


2 2 
(Eia Т A E; n) = En + 2 Ἐπ Δ ο. 2-7 
Substituting 3-6 and 3-7 into 3-4 gives 


2% . . 2 2 Ë _ 
ІА + ІДА + Са А + С ДА + 30+.) E (ға) RS S f(A)f'(A) AA 


EE. 2 ‚_ ДЕ Ы 
S(tg) 5 (cs) Ё іп * 2 S(tg) S (cs) ips А іп 3-8 


Subtracting the perturbed equation, 3-8 from 3-3, the 
following linear equation with constant coefficients is obtained for 


every Eig anc its corresponding rotor angle A. 


I ÁA « Ca AA + 2 S(tp) S*(es) £04) £'GO) ДА 22 506) S“ (es) Ein À Ei, 
329 
Rearranging 3-9 into its conventional form and comparing 
the coefficients with those of the standard second order liner differ- 
ential eauation as was done in Chapter II, give the following two 


equations. 


1 a ο 
АА + ДА «AA 
е ООЛ ώς τα. so ТЕС и ” 
25 (+=) S 5- f(A) £'(A) 28 (42) 5 Gea) f(A) f'A) 
I Bin A Ein 3-10 
=. f(A) f'(A) 
1 AA3420D) ДА+ ДАА = Ein Ein 3-11 
м2 Wn f(A) £°(A) 
Thus 5 
W = 2 50 - (св) f(A) £'(A) 3-12 


ала 
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= Pe DER í τ᾿ 
(08) (оузз) = 2 2 S(tg) S^(cs) f(A) f'(A) I -- 


Tne fact that f(A) and f'(A) occur as products will be very 
important in determining tolerances on the dynamic constants of the 


system. 





We shall now check the above analysis for the iceal system 
described in Chapter II. 


Че have shown in Chapter II that 


Ед“ = S“(sr)[e2;e] 5(sg) [Ase | S 3-1 


Непсе 
1 
Eo = S(sr) ЕЗІ 5 (вв) [Ase | ë 3-15 


Thus, in the ideal case, we see that 


f(A) = S(sr)|e2;e| D [5e] А 3-16 
апа 
Ғ'(А) = ec = S(sr)[e35e| 5 (зв) [зе 5 = 


and finally, we have 


f(A) f'(A) 


5 32 г) u (с) [Ase | 3-18 


which, when substituted into 3-12 and 3-13 gives 






Ste) S (es) (sr) St. [ΑΜ әлә 


1 І 







Ul. Ее 
I (S(tg) S*(es) S*(sr) 5(sg)) 


ба 1 3-20 
N Ttr] 


(28) (оузз) = 5” 


" 
№ 


2(ovss) 
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Equation 3-19 and 3-20 are identical with those developed 
in Chapter II for the ideal system using the conventional approach. 
As mentioned before, the W, and ОШ) (буаз) corresponding to the ideal 
system will be referred to as Wh(ref) and (DR) (ovss) (ref)* 

Feturning again to 3-12 and 3-13 and solving for the product 


f(A) £'(A), we have for the two equations 


2 
f(a) f(A) = —— = кум 3-21 


2 S(tg) 9" (cs) 


where 
I 
2 Site) (сз) 


2 


“9 1. 
O 8 5(+р) 5<(сз) 1 (08) (сузв) ^2 (ТЕ) (оузз) d 


where 2 


ба 
27 Β5{ερ) σος) 1 

We shall now have to assume arbitrary tolerances on the 
undamped angular frequency and on the damping ratio of the system. 
Having made these assumptions, it will now be shown that certain 
definite limits can be set on the non-linear device used in the feed- 
Баск loop with respect to its departure from an ideal square root 
generator. We shall first discuss the Wy, parameter and attempt to 
determine the characteristics of the non-linear device from tolerances 
оп Wg. 

Let W, of the non~ideal system vary from Wh(rer) by some 
fraction of Wa(ref) such as £ f W.(ref)* Thus we are allowed two 


limiting angular frequencies defined as 
Wn(L) * "n(ref) - f Wn(ref) 3-23 
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the lower limiting natural angular frequency апа 


Wn(u) = "n(ref) * t (ге?) 3-24 
the upper limiting natural angular frequency where f is to be defined 
as the tolerance to be set upon the system expressed as a fraction of 
Wn(ref) > 

Thus, we see that for a given tolerance, the right sice of 
3-21 15 uniquely determined for the two boundary values defined in 


3-23 and 3-24. Solving 3-21 for the slope of the (A, Eo) curve, we 


ш f(A) BI S(tg) Sen) f(A) ч: 


and hence two limiting slopes are determined at every point in 9 
plane in which A is the abscissa and f(A) = Ес is the ordinate. We 


define the lower slope at any point (Aj, f(A,)) аз 


`2 
_ _I (Wnlref) - f "n(ref)) E 
E TES fü) Ñ 


2 Ῥ(ύς) μ22Μ] S* (es) le; | 


and the upper slope at any point (Aj, f(41)) as 





они = бей t f (тег 1 иЗ 
2 S(tg) [425M | S (es) [esi] f(A) 

Fig. 3-3 shows the mapping of the A, f(A) plane for a given 
tolerance of W, . Many significant points of interest can be gained 
from a close study of Fig. 3-3. It is first to be noted that the upper 
and lower limiting slopes along a line f(A) = const. are constart. If 
the tolerances on W, are increased, that is we are allowed to deviate 
by a greater percentage from Wn( ref) throughout the operating range, 
the differences in slopes at any point are increased and hence the 


openings between the lineal elements are increased. Such a diagram as 
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POSSIBLE CHARACTERISTICS OF MON-LINEAR DEVICE 
IN ORDER TO STAY WITHIN ASSIGNED DYNAMIC TOLERANCES 





shown in Fig. 3-3 is celled a linesl element diagram? Paths (1) and 
(2) of Fig. 3-3 follow the steepest slope at any point and hence are 
solutions to the differential equation 3-27. Path (3) on the other hand 
follows the lowest slope at every point in the A, f(A) plene starting 
from the same point as (2). Path (3) is a solution of 3-26. There are 
an infinite number of solutions to each equation, each solution differing 
by a constant. Path (4) shows how a path may vary and still give an W, 
within the given tolerances of the system. Thus, it can be seen that a 
non~linear device which has an input-output characteristic whose slope 
at any point stays within the specified tolerances of Wp will be a sat- 
isfectory device as far as the dynamics on Wp are concerned. Те 
(DR) (ovss) tolerances will be discussed in the next section. 

As f(A)—»o, f'(A) OO and of course it will be impossible 
to obtain a device that will start off with an infinite slope. Thus, 
for low input signals, the dynamic qualities of the system will have to 
be relaxed since the product f(A) f'(A) cannot be mainteined within the 
required limits. That is W, will be decreased and the (DR)(ovss) will, 
of necessity, increase. The best we can do is try to obtain as steep 
a slope as possible for low inputs and specify the input below which the 
dynamic tolerances are exceeded. 

We have shown that for the ideal system 

1 m 
f(A) = Eo = S(sr) [e2;e] S2(sg)|A;e| 5 3-28 

which is the equation of a parabola. Let us re-define the equation of 


the parabola as 


м 


f(A) = Eo = PA 3-29 
and denote P(ref) as the constant associated with the ideal system. 


The limits on P are to be determined from the given tolerances on Wn 
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ог (DR) (vss) 
t'(A) » ipa) 3-30 


and 


N 


f(A) (А) = > 3-31 


Since we have shown that 


2 5(іг) 5% сау 2(А) 2" (А) 
Wn = - 


In this case 


| ‚2 
5 5 
Wn = x S(tg) ° (cs) P 2-32 


1 


a 1 
(DR) ( оузз) 215 | 2 5(1ϱ) 5“(са) 1 ҒА) £'(A 


Са | 1 “4 
(DR) (ovss) = > δεις) ες) Í P 3-33 


Taking differentials in each case, 


d Wn = d P 3-34 
and, dividing 3-34 by 3-32 gives 


dW - dP 3-35 
Wn Р 


леи d P E 
d (DR) (ovss) = 2 т S2 ss) Τὶ, ΒΟ 3-36 


and finally, 





also 


а (DR) (ovss) - -ЧР 3-37 
DR) (Ovis Р 
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From the preceding equétions, we have deduced the important 
fact that a fractional change * f Wa(ref) equals the same fractional 
change in the constant associated with the ideal square root curve while 
the same fractional change in Эр abbreviated * f(DR), produces 
an equal and opposite fractional change in the constant associated with 
the ideal square root curve. 

Fig. 3-4 summarizes the results of the previous paragraph. 


Т 


Let fj be the specified tolerance on Wn(ref), i.e. Wy can have any 


value between Wn(ref) гү Wn( ref) апд М - fi Wn(ref)' Correspondingly, 
the constant P will ШЕШЕ any value between P(ref) * ІР (геғ) апа P( ref) c 
fiP(ref). In addition, the COR Безо) can assume any value between 

(DR) (ovss) (ref) in (DR) (ovss) (ref) and (DR) ( oyss)(rer)* fy (DR) (ovss) (ref): 
It is important to note here that for a certain given tolerance on the 
natural angular frequency, there automatically corresponds the same 
tolerance on the damping ratio parameter through the upper and lower 
bounds of the square root curves. 

In Fig. 3-4, we have shown how tolerances on the dynamic 
parameters allow an upper and lower bound of square root curves about 
the ideal curve, Eo 7 P(ref) AR, It remains to connect the development 
in connection with Fig. 3-4 with that of the lineal element disgram of 


Fig. 3-3. We have already shown that 


Eo = f(A) = РА 3-38 
f'(A) - 5 д 
. LP? 
7 2р А5 
2 
E Кер m B 
5 210) 3-39 
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LIMITING SQUARE ROOT CURVES FOR A GIVEN TOLERANCE 
ON THE DYNAMICS OF THE SYSTEN 
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Where we have defined P as the upper and lower bounds of 
the square root curves corresponding to the tolerances imposed upon the 


dynamics of the system. That is, P will assume the two limiting values, 


Pi 7 P(rer) * fi * (rer) 2-50 


P2 * P(rer) - f1 P(rer) E 


Thus, since 


= 1 ae 
P (ref) F "(sr) |. 2-42 
1 1 
Р - 5(.) 5(-4)2 2 4 (у)? g 
2 η 1 : 
diis S(sg)* 5 11 (sr) "(sg)* 37h 
Using 3-19 
2 f) W ° 1 3-1,5 
РИ ry 2 ty р Po EET RERO TM 5 
В о 
Hence, going back to 3-39, we obtain 
Я І 1 3-46 
ЕЕЕ > біег) 2 Ғл hilei) ης 
2 S(tg) Ὁ (56) f(A) 
But, by equation 3-27, we have already shown that 
I 2 
ΠΛ - Wn( ref) 2f Wn(ref) 3—7 


2 S(tg) Э(св)^ 


From the fact that 3-46 and 3-27 are identical, we conclude 
that the tolerance limiting square root curves shown in Fig. 3-2 
provide the limiting slopes not only within the band of square root 
curves but automatically pive the slopes necessary to mep the entire 
plane as shown in Fig. 3-3. Another important consequence in the 


development of 3-46 15 that an assigned tolerance on the natural 
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undamped frequency, simultaneously fixes the same tolerance on the 


damping ratio. 
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CHAPTER IV 


Transient Analysis of the Orthogonal Vector Summing System 


In this chapter, we shall show in detail the response of the 
system to input step and ramp functions. We shall take as our starting 


point equation 2-18, and let E;, be unit step function 


Thus, 
Wn^ até Wn dt 


where we have let E^in be the square of a unit step function applied at 





t = 0. Assume the system at rest at t = O. 
| The response equations to be analyzed in this chapter will be 
solved by the Laplace transform:tion method. The Laplace transform of 


a function f(t) will be expressed as 


о 1200) = F(p) = f «oet =? 


о 
Taking the Laplace transform of both sides of 4-1, and 


assuming zero initial conditions we obtain 
p JE 2(DR) OVSS А | Е 2( ) at НВ 
іш р н р о ‘Р’Пр 
Here, we have designated E^ (p) as the Laplace transform of 
2 
E^ (t). 
Rearranging 4-3, 


1 =, 


2 =y? 
E — N а t 
ote) D  »(p? « 2(DR)(ovss) Wnp * πο 


There are three cases of equation 4-4 to be considered. The 


three cases to be considered are for (DR)(gyss) <1; (DR) (ovss) = 1, and 
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апа (ГЕ) (оузз) >1. 

Before proceding further it will ve important to note that 
because we are dealing with a second order systein with a finite moment 
of inertia, it is impossible by any finite forcing function applied at 
t = QO to produce a displacement of the shaft which will in turn produce 
an Eg(t). This is of course true even though we hzve made the simpli- 
fying assumption that there are no dynamics involved between the shaft 
and Eg. Lot only will there be no displacement at t = O to a forcing 
step function, but the rate of change of response is also zero. Thus, 
the square root of the time response to te derived from 4-4 will nave 
to setisfy these conditions. 

Equation 4-4 when solved by the usual method of partial 


fractiors gives as the time response to a unit step function o? Big, 


- : N 
τ) = ] — cep ον е πον.» nt ring 2-920... 


V je λαο 


~(DRA) (oyss)) (PR) (ovss) <1 4-5 


where the dynamic response angle, abbreviated 


` 2 
Е l - (DR 
(DRA)( ` = tan 1 V 17. (R^ (ovs) 
0У55/ (DR) 
(ovss) 
To find the initial values of ESSE) and its first three 

derivatives, we note the following from "Transients in Linear Systems", 
by Gardner and Barnes." 


Given а {о} - Е(р) 


Тһеп 








f(O+) = lim p F(p) 
pot se 
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= f(O+) = lim Е F(p) - £(0+)p| 


P oo 


42 


5 £(0+) = lim E F(p) - £(0«) р? ~ E г(о+)р/ 





pl» 5Ὁ 


Э а 42 N 
Ф” f(O+) = li ERE τη..." 5, 
TE (0+) im | pt F(p) - #(0+)р m (0+)р 2 /р 


| р-—— °° 
Hence, applying the above criteria to E^. (p) developed in 
4-4, we obtain 


Εο΄ (0) - 0 


2 
S. ES(0) s Ww? 4-6 
dt? 


ο ο 2 
Е Ео“(0) = - ОЕТ) AR 


We shall abbreviate (DR) SS) to (DR) for the remainder of 
this chapter to shorten the notation. 


From 4-5, we see that 


ose ede 1 е(09 96 аяп(ы, _ (рн)? % - (ВА) (оузз)) 
V1 - (DR)* T 


Equation 4-7 represents the response of a passive system. 


$ 


Thus, even with no damping in the system, the maximum value of the 
second term on the right would be unity. Since a fairly large amount 


of damping will be used for the system, the second term of 4-5 will 
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rapidly become small with respect to unity, and thus enable us to 
approximate E9(t) closely by two or three terms of a binomial expansion 
of 4-7. 

It is shown in "Differential and Integral Calculus", by 


Granville, Smith, and Longley that the expression 


m(m - 1) х2 + . 4-8 


(ых) = Deseret qo e ө ο 
and is valid when m is not an integer providing |x| <1. 


Applying 4-8 to 4-7 we obtain as a good approximation of 


Eo(t), 


E(t) = 1- 1 en“ (DR) Wat sin(w ,/1 = (DR)® t - (DRA) 
2 [1 - (08) т (ovss)) 


-2(DR) Wht 
ΕΞ G- NE sin Gi y - (DR* t - (DRA) (ous) 


4-9 
Taking next the case where (DR) = 1, 4-5 becomes 
2 2 Į 
τι ο το SS 
o «P n р(р 5 Wn)? 
Taking inverse transforms of 4-10, we obtain 
E 2(t) = 1 - e""n* (Wit + 1) 
ος, finally, 
Eo(t) = h - eT n* (ut. | 4-11 


Expanding аз before, in a binomial series, 


Eo(t) = 1 - jent (Wut + 1) - je ^ n^ (ut + 1)2 - . . 4-22 


Taking finally, the case where critical damping is exceeded, 


we have (DR) > 1 and hence the characteristic equation associated with 
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4-5 factors into 


(p * (DEW, - Wig (OR) - 3p (м + a, (0R)? = 


Thus, 


и < 
р(р + (DR)W, - W, (DR)? - l)(p + (DR)Wn + “п ү/(юв)* жү 


Тһе time expression of 4-13 becomes 


r. | DLP 
Eo^(t) = 1 + š Cm? an JE -ici ο. + Wr (DR) EIE 
- (DR DR)“ = 


1 .- aw, - w, DR)? - 1 { 


MESCUEM p 
(DR)? + ш - 1-1 и 


Eo“ (p) = 4-13 


~ 


+ 


Thus we have finally, for (DR) > 1 


(К = Ё +} PT es СЕЕ Қою, + Wr JR)? - || L 
(DR s (pa pu Ἡ 


ЖЕТІ == „о те iy fon)? = 2 fe š 
πο... 


(DR) > 1 4-15 

It is not feasible to expand this expression by the binomial 
approximation. 

Curve I of Fig. 4-1 shows the solution of 4-5 for (DR) = .1. 
If the square root of tne values representing curve I are taken at 
every point, we have the exact solution shown as curve II, and which 
represents 4-7. Since this latter result is difficult to obtain 
analytically, we have expanded 4-7 through the binomial series approxi- 
mation to three terms in 4-9. Plot III shows the approximate value of 


Eg(t) to a unit step input function using only the first two terms of 
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4-9. it is to be noted that the approximation is very good except 

for the initial period of T = .1. Since damping ratios higher than 
.1 are actually to be used, the approximation will be even closer to 
the actual solution except that as the damping ratio increases the 
time for the initial overshoot will correspondingly lengthen and hence 
the initial large deviation will last for a longer value of - 

Throughout this thesis, it has been emphasized that the 
terminology used such &s damping ratio, and natural frequency refers 
only to the linear second order system relating the squares of tne 
output and input functions. However, one must be careful and not 
apply this terminology to the specific response &ssociated with the 
input and output functions themselves. 

In this presentation, we shall use the πη. decrement 
as our basic description parameter for the system for (DR) < 1. Loga- 
rithmic decrement will be defined as the natural logarithm of the ratio 
of the amplitudes of two values of the response to a step function 
differing by a period as for example two successive maxima of oscil- 
lation as shown in Fig. 4-1. 

Referring to Fig. 4-1, we let A, and Aj be the amplitudes of 
two successive peaks of the linear second order response of the squared 
functions while A,' and д! are two successive peaks of the response 
for signal input as a function of signal output. 

Hence, the logarithmic decrement for the linear second order 
response, abbreviated (LD)(jin) * log =, 4-16 

From Equation 4-5, we see that the instantaneous height of 


the overshoots shown of Fig. 4-1 is 


1 _ | DIE c A 
т тй” e- (DR) Wnt ONE - (DR)? t - (ПВА) (оузз) 
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and hence using 4-16, 


(12) (1: = log Т Ж ү (DR) Wp T е1 
їп) 7 (DR) Wyte 


since A, and Ag were taken at the peaks and T = to - t,, the period of 
the damped oscillation. 
From 4-17, we obtain 


(LD) (315) = (DR)W, T where T is the time for two positive 


2м = 277Г 
ссезз1уе peaks. But T= = “° ———P"O thus 
Su p М ИЕ (DR) 9 


+ 2 7 (DR) 1-18 


(12) 
(lin) E _ (DR)^ 


We shall next defire the logarithmic decrement for the actual 


response of the system to 8 step input as, 


= Al = Al Ей || 
τυ. с V A> 2 s 54-19 


From 4-19, we nave immediately that 


(LD) (lin) = 2(LD)(ovss) 4-20 


Using 4-20 in 4-17, we obtéin the important result that 


ο a ee η 4-21 


ТШ ІНЕ 


This last result, gives a neans of describing the non-linear 
response in terms of a parameter describing the linear squared system. 

Up to now, the development of this chapter has been with the 
ideal system which assumes a perfect square root function in the closed 
loop. The last chapter showed how close a square root function had to 


be generated in the closed loop in order to stay within certain pre- 
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assigned tolerances in the damping ratio and the undamped natural 
frequency. Returning to this line of argument, if we have a change 
іп the damping ratio, d(DR)(oyg5) then from 4-20 we obtain 


_ Т(1» (рв)2) E 
d(LD) (oss) Do Г d(DR) (DR)< 1 4-22 


Applications of the response of the system to a forcing step 
function will be made in connection with the laboratory test system. 
The constancy of the dynamics for various magnitudes of the input signal 
will be tested by applying small step function increments upon the dif- 
ferent static operating levels and recording the output response. If 
we have a perfect square root generator in the feedback loop, the 
response to a step function increment will not depend upon the magnitude 
of the forcing function, and hence a non-dimensional response curve can 
be shown by plotting Bo for any point of the operating region.  How- 


Ein 
ever, we shall use an imperfect square root generator in the feedback 





loop and hence small increment forcing step functions will be used and 
the operating points specified. Chapter VI will show in non-dimensional 
form, such responses as obtained from a Sanborn recorder. 

The previous analysis with πο ο decrements was made on 
the basis that the (ТЕ) ( әузз) <1. We have not developed an easy рага- 
meter to describe the non-linear system in terms of the linear squared 
system for (DR) (ovss) Ξ 1. However, the response can be obtained 
fairly easily from 4-12 or 4-15. 

It was felt that to make this chapter on transient analysis 
complete, an analysis of the system response to a ramp forcing function 
should be included. Much study was devoted to the comparison of the 


response of a regular second order system to the one presented here, but 
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the results were disappointing because no simple relation or design 
parameter could be found relating the two. However, a brief develop- 
ment of ramp responses follows for (DR) < 1. 

Taking the Laplace transformation of 2-18, and assuming zero 
initial conditions, we obtain 


2 2 
ре а a MEE 
n n 


or 


2 
Ein (p) 4-23 


a p^ + 2(DR) р + 1 
Яр Wa 


Eo (p) = 


In developing the response to ramp functions, it will be 
best to work from the point of view of errors (both steady state and 
transient) between the input and output functions. 

Treating 4-23 first as a linear second order system in DM 


we define the error for such a system as C m ως = E^ (р) 


where 
1 2 2( DR) 
2 5; Ein(p) Luge Р? + πο P| 
Eo^(p) - E*in(p) - zn 
o «P in\P 
a p^ 4 2(DR) р+1 
Wn и 
ог 
P». «οὔ E2;n (p) p? * 208, p]. 
Eo (p) - E in(P) 4-25 


p? + 2(DR)W. p + Wx“ 


τι ο ο a unit ramp function = t, 


И 
Е? (р) Б 
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 —F sn u μμ LL κ 


апа the linear second order error equation becomes for DR < 1 


p + 2(DR)W, 
Eo (p) - Ein (p) = 





p(p* + 2(DR)W, p + W,“) 
Ν p + 2(DR)W 
p(p + (DR)w, - jug, /1 - (DR)2) (p « (DR)Ws « jug /1 - (DR)^) 


4-26 


Ко К К 
= ———À рае 
P р+ (08) - М1 - (08) р + (DR)Wy + 4/2 - (00 


where 
Ky = (р) (Е = EE) т алы 
п 
р= 0 
Ир < (pR) = Wa. Ho (DR)“) (E 2 - P) 


р = ΠΠ 434,4/1 = (DR)^ 
"entem doe 


] tan 
ο. 


2 jW,./1 - (DR)? 


Hence, the complete time expression for the linear second order 





error response to a unit ramp is 


Е,2(%) - Ειρ΄ (ο) Ξ AM + AG ет (08)! nt sin (Wa A _ (ов) а 
(DRA) eves) 5-21 


where 





(DRA) (oyss) = барта 2 Vi - (DR). (DR) « 1 


(DR) 
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If, on the other hand, we let E;, = a unit ramp function = t, 


then 


Bin’ (P) a М z 


and 4-25 becomes 


2) _ в. 2) = р 20 WEB, 


P^ [p + 2(DR)Wpp + is? | 


2|p + 2(DR)w, | 
рр * (Юм, - jW,,/1 - (DR)? || p + (DR)W, + Sif ~ (DR)? | 


(DR) X 1 eae 





Expanding by the same method of partial fractions as in the 


previous case gives 


Е 211 - pa)?] , 4(DR) , , 2 
ΤΠ : МОЛ мул. πα 
4,2 Wh Wa~(2(DR)2 = 1/2 - (08)? 


ον sin(W,,/l - (Dg)? t + (DRA) ұза) ΡΕ «1 





4-29 
where 
\/ D 
(DRA) a ND] хов) 1 = oe. ms 
(оуз5) (DR) 2 


An interesting result from 4-29 is the fact that if the 
Da г is made equal to .5, then the first steady state term in 
4-29 is eliminated. However, the second term shows that the error 
between the scuares of the input and output builds up linearly with 


time. 


60 





CHAPTER V 


Steady State Harmonic Analysis of the 
Orthogonal Vector Summing System 


An important method used to determine the dynamic charac- 
teristics of a general electro—mechanical system is the method of 
steady state sinusoidal harmonic analysis. Amplitude and phase graphs 
are plotted as functions of forcing sine wave frequencies. In this 
method, all transient responses are considered as having died out 
leaving only the steady state conditions. The dynamic quality of the 
system is then determined by the pre-assigned limits set upon the 
magnitude change and the phase change in comparison with some reference 
magnitude and phase determined at some reference frequency, the latter 
usually being zero. 

In the system being considered here, it is impossible to 
follow such a procedure because of the squaring properties of the 
torque generators. Thus, only the absolute value of a function acts 
on the system, and of course (соз wt |, which looks like the following 


function, 


Wt —- 


has too abrupt & change to be of any value in evaluating the speed of 
response of the system. 
If, on the other hand, we impress upon a constant input F, 


a function of Wt, F4 cos Wt, where F4 = Ес» the discontinuity can be 
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prevented. Thus, let Ey, = Fo + F] cos Wt. Since the equation giving 
E in terms of S is linear, all of the solutions for inputs of Eip 
will be found through the squared functions. The equation for Eo? in 


terms of Ein is 


2 Р 
і а 2 2(DR) а 2 ο. - 2 
— — E — E“ +E = Ejn 5-1 
Wn~ dt c dt 


An important significance of this linearity is that if Eip 
consists of a sum of terms, a solution of the equation can be obtained 
for each of the terms independently and the overall solution obtained by 
superposition of the separate solutions. 


Taking an input function, Ey, = Fo + F) cos Wt, and squaring 


gives 
Ein" ж Ес * 2 Fg F1 cos Wt + FA cos “Wt 5-2 
From trigonometry, we have the following relationship: 
cos^Wt = $ + À cos 2 Wt Pr 
Th 
из, το 52 
ah Кт о le 
Ein = Fo + 2 Fo ἡ] соз № + 5 2 cos 2 Wt 
or finally, 
2 Е- 2 
Е 1 
2- 2 — -- 
Ein = Р + 2 ΤΣ Fy cos Wt + 2 cos 2 Wt 5-4, 
F^ 
Taking the p * 5 . erm and inserting it into equation 5-1 
p 
gives the response ies = Fo + T in the steady state. It is to be 


stressed that this harmonic analysis deals only with steady state 
responses of the system. 
The steady state solution of equation 5-1 to the fundamental 


sinusoidal forcing function will now be shown. 
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Rewriting equation 5-1 as follows: 





2 
Е О ο cen е5 
Wn dt Wn dt 


Re 2 Fy Fy 5-6 


Let 


S 


Assuming & particular solution of the form, 


- JWt 
2р А е 
2р = JWA eit 
2р = - We A еле 5-7 


Upon substituting the above functions into equation 5-6, 


we obtain 


NE OR) ο o pop; 
n Wa 


ҚЫ 


T 


ОКЕ 
о 


ο 533 
T J ( D 


n n 
Thus 


2 2 ЕС Bi e JWt 
1- 2 + j 2(DR) W 


Wn Wn 


ace Re 2 Fo Fj (cos Wt + 3 зіп Wt) 5-9 
1 - (FR)* + 3 2(DR)(FR) 


where 


|= 


(ЕВ) =; 


т, 
> 


63 


_ ве 2 Fo Fi(cos № + j sin Wt)(1 - (FR)* - j 2(DR)(FR)) 
ΠΝ ЖШС CNN 2 
Е z (кв)? | + [2(0ς) (58) | 


ог 
| 2 Fo Fy (1 - (FR)*)cos Wt + 2(DR)(FR) sin Wt 
а πας 
(1 = (FR)? | + [ 2(DR) (FR). 
5-10 


or finally 


2. 2 Fo F} 
yp- 2 ~ роз |а 


(ЮВА) у = қап - 7. 


On a similar manner, the response of equation 5-1 to a 


2 
function c cos 2 wt is of the form 


cos (Wt « (DRA);) 5-11 


where 


2 
ИНЬ (Wyt + (DRA)>2) 5-14 


vh - (83? |? е |2008) (вв) |2 


мћеге 


san1 - 2008) (Е8)] 


с 1- (ғн)12 


Since (FR)) = = = 2 - = 2(FR) 
п 


Thus, 5-12 becomes 


2 
EET е μα cos (2Wt + (DRA)>) 5-13 


Vl. - в)? ]? „ [(ов)(тв)] 2 
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where 


(DRA). - tan. E: (Ән) (ЕН) 


ШЕТІ ЕБ 
In the case of 
L(FR)? «« 1 
(DRA), = 2(DRA); 5-1, 


This can be seen from the fact that for (FR) << 1 the tangent 
of the angle can be replaced by the angle and (DRA) > = - L(DR)(FR) as 
compared to (DRA);* = - 2(DR)(FR). 

This approximation is good even for (FR) = .3. To show this 


fact, we shall let (DR) (ovss) = 


20 
l =- .09 


(DRA) 1 tan71 - 


RE 33 S 


Letting (DRA), = 2(DRA),, gives - 66.8? for (DRA) 5 


Actually, 
(DRA), = tant E = ре 
Thus, corresponding to 
Ein? = Fo* 5 + 2 Fo Fy cos Wt к cos Wt 5-15 
we have in the steady state 
Ед? = Rae (ee Е qia cos (Wt + (ОКА) у) 5-16 


"M (ға)2|2 + [2(рв)(ея)| 2 


2 


н ————— cos (2Wt « (DRA)2) 


vi 2 A(R)? | 2. |, (ot) em? 
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where (DRA), and (DRA)2 have been defined in 5-11 and 5-13. 
It is now necessary to relate in some fashion, the squares 

of the inputs and outputs to the actual inputs and outputs themselves. 
Starting with the expression Eg? - Ejn? - (Eg + Ein) (Eo - Ein), 


we can make the approximation 


Eo? - Ein? - 2 Eig(Éo - Ein) if Eo - Ein. 5-17 


This approximation will be valid because we shall restrict 
ourselves to small deviations of Eg. 


Dividing both sides of 5-15 by ps gives the result 





2 
Ео =e GS 52 = - || 5-18 
Ein’ Ein 


we shall define the deviation of the performance ratio of 


the orthogonal vector summing system as 


2 


Eo = 1 s - 1 
DOR) CSS) E 2 Е, 12 5-19 
2 
where 5-15 and 5-16 are to be used to give the expression for T. 20 l. 
in 


Information from 5-19 will be broken into two parts, each 
part conforming to certain specifications. Іп this particular treatment 
we shall use the damping ratio associated with the linear second order 
performance equation relating the squares of the input and output 
functions as the necessary parameter to give us a family of performance 
ratio curves plotted as a function of frequency. The first of these 
aforementioned parts will be a comparison of the ratio of the absolute 


2 


magnitudes of Eo? and E;,~ while the second part will give the phase 
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| 
| 
| 
| 


2 2 
lag of Eg as compared to Ein as a functional forcing frequency. It 
is to be noted that the amplitude deviation will be a comparison of 
the maximum amplitudes of input and output quantities. From the infor- 


2 


mation obtained from the ratio of the absolute magnitudes of E,“ and 


Din 


we can then use 5-19 to find OCE 5) as a function of frequency. 
In the expression Гог Еш, that is Ey, = Fo + Fy cos Wt, Fo 

апа Еу һауе been specified as not necessarily equal. This general 

approach has been carried through with an eye on the future applications 

of this type of analysis, In considering the quality of response of 

the system, Fy will be made equal to Fi. However, in considering the 

effect of noise on the system, F4 will usually be a small fraction of 


Fo. We shall apply our analysis only to measure the quality of response 


of the system, considering a study of noise beyond the scope of this 


thesis. 
Ihus, letting 
Ein = Fo + Fo cos Wt 5-20 
Треп 
Е. 
“in = 1 + cos Wt 5-21 
Го 


Let Fy be redefined as Ein(ref) 


Pin ντ созе 5-22 
Esn(ref) 
2 
|I + cos Wt 
Ein | 
E in (ref) 
О TT 217 
Wt —» 
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Тпеп, 


2 
Επ. = 1.5 4 2 cos Wt « .5 cos 2Wt 5-23 


Ein^(ref) 


After substituting this quantity into equation 5-1, 


Eo SPARE — sk enn er = cos(Wt + (ОКА) ч) 
Е 
in(ref) vh = (FR)2|2 + [2 (28) (ονος) (FR) | 2 

3 ‚2 — R cos(2Wt + (ΑΛ) 2) 


vh 4 (FR)?| Я ао) ος N 


where, for the sake of completeness 


_1  2(DR)(FR) 


tan - 


(DRA)1 


1 = (FR)- 
E DR)(FR 
1 - 4(FR) 
E? 
Fig. 5-l shows how the input function 00 —— of 5-23 is 
E2; E in(ref) 
compounded. It is to be noted that has a frequency, W, but 
in(ref) 
the required addition of the second harmonic causes distortion as shown 
Е,2 
in Fig. 5-1. For sufficiently high forcing frequency, м 
2 Е2. 
5 іп(ге”) 
lag the input as shown and the absolute value of —5— ———- will drop 
in(ref) 


off in general or it could very well peak for low damping ratios near 
the natural frequency of the system. 
Fig. 5-2 shows in non-dimensional form a plot of the 
ШЕН) (2/55) as a function of forcing frequency using equation 5-19. 


E2 
The absolute value of —*2——— can be seen from equation 5-23 to be 4. 
in(ref) 


68 





Ф| 
о) 
| 





ΠΕΠ 
ИМ з0э2 7 \ u 
CRT CCE 
YN Е ΜΝ T [| 
"ЯСЫ 
ο. ο EE NAN 


τει ДАТА) 
ΤΗ ΣΤΗΝ ΤΗ ΤΗΝ, 
_ у үү 


пейш z 
^ 
2 


< 





λα AJ (зал) 3 d али š LA » 
ШЕ N V 
01 лз 30 NOSIJVdNO9 (NV 1ñ dB! в" NV 
(13110 (3711ΥΝΝΟΝ 3Η: 30 №01 11509009 
amd AÑ 





— OILVY SONLIIGNYV 


69 


ΟὋΙ 


ШШ N ШЕ [ 
yd) = [= 
TT 









ΠΠ ΠΠ T 
ШІП ΠΠ 1η 
А И T 
НС 
ΠΠ ЫЕ ТТТ 
ТО ЕИ, 
ШШ ШШ 












ТА | wo 
JNINNS 0123) em 
ШИ MALI 30 35009535 03711УМ0І5НЭМІ0-КОН 






0-5 380913 


ШШ, КТЫ ШШ | 






A LL ΠΠ η lu d 












Ша 
(492) "3/03 


m ча o » 





70 





2 
Е 
ο 
2 can be calculated from 5-24 with (DR) (oy ss) as the parameter. 
in(ref) 
The second harmonic component will cause a resonant point at (FR) -.5 


in addition to the usual resonant point, hence the broadening of this 
peak over the usual form in which it is seen. The fact that for low 
damping values the resonant peaks become narrow and high, brings out 
an interesting point in connection with the curve for πη к 
Неге ме have two distinct peaks, the first due to the second harmonic, 
coming at a frequency ratio slightly less than .5. 

In order to find the phase response of the system to the 

input ince on αι, = l + cos Wt, it will be best to look at the 
Bin(ref) 

squares of the outputs and inputs in rotating vector form. In setting 

up the relationship between the input and output rotating vectors, let 

us assume (DRA); = 30° and (DRA)> = 60° as a typical case. Further, 

let us assume for simplicity that the amplitudes of the response vectors 

are the same as the corresponding input quantities. 

It is obvious that the phase angle between input and output 
will not remain constant throughout a complete rotation of the funda- 
mental vector due to the second harmonic vector rotating twice as fast. 
This changing of phase angle throughout one period becomes more pro- 
nounced аз (DRA) builds up. The best that we can do for the given 
system is to predict where the (ΒΑ) (ονςς) starts to build up to a few 
degrees. Inasmuch as the fundamental vector is four times the second 
harmonic for most cases (not near resonance for low damping ratios), we 
can, for small phase angles, approximate the total phase from the funda- 
mental as can be seen from part I of Fig. 5-3. We also must define this 
phase angle to be determined at t = 0, 2 , εἴς. 45 shown in part I of 


Fig, 5-3. Fig. 5-4 gives the phase lag of the fundemental as a function 
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 ——  - GONSTANT INPUT—OUTPUT ASSUME: 





HeH = INPUT (DRA), = 30° 
= OUTPUT (DRA), = 60° 
------------- =o- 






( DRA} ovs S) 





FIGURE 5-3 
INSTANTANEOUS RELATIVE POSITIONS OF OUTPUT AND INPUT FUNCTIONS 


FOR THE FORCING FUNCTION, <--=!5+2coswt+Scos2wt 
intref) 
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of frequency. Caution mist be taken to use Fig. 5-4 only for small 
angles. 

In concluding this chapter on the steady state harmonic 
analysis, it must be pointed out that approximations have to be made 
due to the non-linearities of the system. These approximations give 
good results at the point where the magnitude of the (D)(PR)(ovss) 
just start to build up and where the phase lags are only a few degrees. 
We have assumed in order to use the material presented here that а 
non-linear device has been inserted into the feedback loop so as to 
keep the dynamics essentially constant. 

This chapter emphasizes the difficulties in obtaining a 
method of analyzing a non-linear system so = to give a clear cut 
picture over the entire operating range. 

It is felt by the author that a step function response type 
of analysis is the more suitable method of ascertaining the speed of 
response for this type of system. The experimental results presented 
in the next chapter will bear out this opinion. This chapter was 
included mainly in keeping with the underlying spirit of this thesis, 
namely that it be a study of the dynamics associated with the given 


systen. 
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CHAPTER VI 


Design and Results of a Laboratory Test System 


The general theory developed in Chapters JI, III, and IV 
will now be applied to the design and testing of an actual system. 
Because of the fact that exact specifications of microsyn units are 
still considered restricted, we shall not use the actual values, but 
we have chosen the specifications so that the test results do actually 
give a true picture of the operation of the system. 

Thus, using Fig. 2-2 as our starting point, we shall use the 
following design parameters for the system to be designed in this 


chapter. 


2 dyne cm 
S(tg)[i2;«]* 42-4 TT 6-1 


volt 


S E - 3.25 x 107 --------------- 
(sg) |A;i,n,e| r&d ma cycle/sec 


і (ех) = 10Ọ ma 
N(ex) = 400 cycles/sec 
A( rot) (max) = 109 = „17% гаа1апз 


I(tsm) = 330 gn cm? 


Using the above excitation values for the signal generator, we obtain 


Е volts 
3 (sg) |е | к гад ns 


Since 50 ma is a desirable upper limit of driving current 
for each torque generator, the maximum torque produced per torque 


generator is 


το 


‚2 
M(tg)(max) = S(tg) [i2;M] i (max) 


hewn бше оп х 2500 ms 


та 


106,000 вупе см 6-3 


Hence, the total maximum torque produced by tle maxium input 
Signals will be 212,000 dyne cm. Thus, the maximum torque produced by 
the feedback voltage Eg will have to be 212,0CO0 дупе cm. 

Fig. 6-1 shows in schematic form the torgue generator winding 
arrangement. i, and iy are connected in single axis arrangement on 
separate microsyn units. ig, the current due to the feedback voltage 
Eo, is connected in single axis arrangement but windings of poles l and 
3 of each microsyn unit are connected in series so that the feedback 
current igo produces twice the torque that it would if it were to actuete 
only one microsyn unit. This is necessary because we need to keep the 
maximum torque generator current at 50 ma. 

referring again to Fig. 2-2, it can be seen that in the steady 
state, 


M) + Mọ = Μα 
2 а 3 Е К 2 
*(es)s[eii] 9681 &  * © (cs)gfesi] “(tala “y => (es) fesi] (te), "Ὁ 
6-4 
Since 3(сз)1 = 3(сз)2 


апа 3(tg)1 = 9(tg)2 


6-5 


T 


2 M cM M ТЫ Е < 
" (es)1,2 "(tg)1,2 5, y | = а); “(ta)3 9 


Непсе 
к. Tf 
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"ur ІІ э * 
"Au uu malli 


T7 


ЭА еды [е] а,ә = (еөз [гл] шю; 6% 


For the system described here, we can define а fictitious 
third torque gererator as shown in Fig. 2-2 which will have twice the 
sensitivity of a single toraue generator. The resistance of this series 
connection of torque generator windings will be double that of single 
microsyn operation. However, because of the fact that a current source 
actuates the torque generators, a change in resistance in the windings 
Coes not change the driving current since current sources have infinite 
internal resistence. 


Referring to 6-6 and noting that 


ο κ ο 
5 es) [esi τα η le;i] 


S(es)3 [e;i] = 707 3(сз)1 αι 6-7 


From equation 2-8 we have defined 


: = | z 2 1 Ж 
S(ovss) [4;X | 7 5(&в) 1] 5 (cs) [e;i] š (sr)|e2;e] "(se) [Ass | _ 
wnere it has been assumed in 2-8 that all the torque generator and 


current source sensitivities are respectively equal. However, in the 


laboratory system we have defined 5(tg)5 |1?;u] = 2 (tg), 2 [2° 3] 
3 9 S 


Nevertheless 2-8 is still valid here since we have kept the product 
S(tg)[i2;x] 5° (сэ) [esi constant. 

It is desired to use the full 10° or .174 radians for rotor 
angle in designing the system, this angle being produced at the maximum 
toraue level. 


Hence 
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- 212,000 dyne cm 
S(ovss) [sw] = 


174 radians 


1,220,000 dyne cm/rad 6-8 


Letting Ex(max), and Еу(шах) be 50 volts each, 


OMS ES 20 amas - 1 ma/volt 6-9 
if Ex(max) = By(max) * 50V 

Eo(max) = x] 2500 + 2500, = 70.7 volts 6-10 
Непсе 

Ῥ(ο5)η = áma (max) = .707 ma/volt | 6-31 


This latter value checks with the result of 6-7. 


Using equation 2-8, we obtain 


СЕН = Lo S(oexss)la;l ____. б=12 


5(tg)5 125] 5° (ев)з [ез1] 5(5в) m 


1,220, 000-29-97 


2 
„58 х 102 Зе сш у ,5 лау у.з и 
та volts rad 


Eo Ес 





2210 volts? 
volt 


or 


volts 
|е, | 7.0 == 6-13 


Figure 3-1 shows the characteristics of the square root generator. 


79 


To actually obtain a device that will have a performance 
characteristic similar to Fig. 3-1 is very difficult. This is especially 
true near the origin where the slope of the curve becomes infinite. 
However, through the use of vacuum or selenium rectifier switching cir- 
cuits in cascade with a linear amplifier, the curve of Fig. 3-1 can be 
approximated to a high degree except at the origin. The latter part of 
the chapter will be devoted to the design of a suitable attenuator which 
will stay within certain specified tolerances of the dynamics of the 
system. 

Fig. 6-2 shows an electrical schematic of the system from 
which the experimental results will be obtained. High gain DC ampli- 
fiers with unity feedback, feeding into parallel connected 6SN7 triodes 
in each case act as current sources. 

Тһе DC amplifier arrangement es shown in Fig. 6-2 can be 
considered current sources because of the following conditions. The 
gain of the DC amplifier is very high with a phase inversion and the 
input impedance is also very high. Thus referring to the figure below, 
the input to the DC amplifier (point G below) is driven very close to 


ground because of the very high negative feedback. 


Ein Eo 





$ t 


R2 


Furthermore, since the input impedance to the amplifier is very high, the 


current through R] equals the current through R2. But since @ is ideally 
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at ground, we can consider the current through Hj as depending only 


on Eig and В] leaving Ro to vary at will. However, through current 


considerations 
Ry Ro 

or 
о 6-1, 
Ein Hy 


Hence, Ej is determined solely by the ratio of . In the 
actual system we have fed Ep into a cathode follower and the voltage 
at the feedback point of the cathode is still determined only by the 
ratio of the feedback resistor to the input resistor. The cathode | 
resistor is then calculated to give a certain desired current for the 
given input voltage. Of course, we have to assume that the output 
tube operates below any non-linearities due to saturation. The cathode Í 
resistors are determined so that maximum currents of 50 ma each are 
delivered for each respective maximum voltage input. 
Since the maximum signal from the signal generator is only 
about 2.25 volts rms (400 cps), an AC amplifier has been provided with 
a gain of around 500. This amplified signal generator voltage is then 
full-wave rectified and filtered before being fed through the non-linear 
device as DC. The output of the non-linear device gives the desired 
solution, Eg. The components within the dashed line box provide the 
5(sr) [е2;е| of the box as shown in Fig. 2-2. Because of the fact that 
the maximum forcing frequencies will be on the order of a few cycles, 
the square root section will be assumed to add no additional dynamics 


to the system. A VR-75 regulator tube is used aroun?’ the feedback DC 
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amplifier in order to protect the amplifier from saturation when testing 
the system in the unstable region described in Chapter II. 

We shall now &pply the development of Chapter III to the actual 
system just described taking into careful consideration the need for the 
stop near the null in order to maintain a stable operating range. 

All of the operating components of the laboratory system have 
been specified with the exception that the square root generator was 


left ideal and specified by the equation, 
I 
Eg = 47.0 El sg)? 


In Chapter IIl, we have denoted the ideal case as 
: i 
Ко = f(A) = P(rer) ^ 6-15 
where 


P(ref) D 5 (вг) [е&;е| S (e) [Ae] 


Numerically 


юе 


Р d volt a volts 
Gen ^ Wo ug n CENE 


ог 


E volts 
ШЕСІ) 527170 ss 6-16 


In setting up a workable system, let us arbitrarily set a 
tolerance on the dynamics so that the natural undamped frequency and 
the damping ratio do not vary from Wy(ref) and (DR) (oyss)(ref) by more 
than 25%. This tolerance may seem somewhat large, but it will be seen 
that the principles underlying the theory developed in this chapter are 
more clearly and easily shown if a fairly large tolerance is assumed. 

The expression for the natural angular frequency of the ideal 


system as given in 2-15 is 
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"n(ref) = | Stoves) ase] 


T 
| ι- "TEL <. Using these 
where S(ovss) [45h | = 1,220,000 dyne cm/rad and І = gm cm °“. Using 
values 
Wn(ref) = 60.8 rad/sec 6-17 
Thus referring to 3-23 and 3-24, we obtain 
Wn(L) = 60.8 rad/sec - (.25)(60.8) rad/sec 
Wn(L) = 45.6 rad/sec 6-18 
and 
Wn(U) = 60.8 rad/sec + (.25)(60.8) rad/sec 
Wn(U) = 76.0 rad/sec 6-19 


The two limiting slope equations are thus obtained from 
3-25 and will be defined as 
Жы 


f(A) 6—20 





ЕТЕЛ wolves € 


Hence, for the case in question, the upper and lower limiting 


slopes become 


8130 
Г! (А = 6-21 
(A), ΤΑ 
22.600 
! = ——— - 
From 6-15 and 6-16, we obtained the ideal square root 
equation as 
lt 
Қо = f(A) = (170 Ж 6-23 


Using the arbitrary 25% tolerances, we obtain the limiting 


square root curves as shown in Fig. 3-3 as 
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volts АЗ 


ο OER VOT a 
or 
1 
RUD SEIS CE À 6-21, 
similarly 
_ volts ,À 
f(A)y = 212.5 тайт А 6-25 


Taking the functional derivatives of 6-24 and 6-25 and forming 
the respective products f(A) f'(A) we obtain 


ї'(А)ү, = B 6-26 


22,600 


(А) 6-27 


f'(A)U 


Equations 6-26 and 6-27 check with 6-21 and 6-22 and hence are valid 
for obtaining the limiting slopes throughout the А, f(A) plane. We 
shall abbreviate A( rot) to A to simplify the notation. 

Fig. 6-3 shows in detail the limiting curves of the square 
root generator with the assumed tolerances of * 25% using the full 10° 
of shaft angle. The limiting slopes are determined at any point A, 
f(A) by f(A), and f(A)y. However instead of fixing the curves at the 
origin, let us fix the upper end point at the ideal values corresponding 
to A(rot) = 10° or .174 radians which gives an ideal E = 70.7 volts. 
If we start at A(max)» Eo(max) and always follow the lowest limiting 
slope at every point, we notice that it is possible to stay within the 
dynamic tolerances and intercept the E, axis at a value of 47 volts. 
This path is of course trivial because it would be impossible to form 
the right triangle addition for any value of the hypotenuse less than 
approximately 47 volts. Starting again at А (тах) › Eo(max) and this 


time following the highest limiting slope at each point we find a much 
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more useful situation. In this case, it is to be noted that it is 
possible to compute Eo down to the limiting value of zero. We observe 

in our case that for A= .064 rad., Ey = 0. The fact that it is possible 
to allow a rotation of the rotor before allowing an output suggests the 
feasibility of using an electrical stop at a small positive angle in 
order to overcome the condition of instability at the origin described 

in Chapter II. 

It is to be emphasized again that the loop gain does not affect 
the steady state conditions of the system except to change the angle of 
the shaft, but does effect the dynamic qualities. Thus, doubling the 
loop gain from that amount shown in Fig. 6-3 reduces the maximum angle 
increases the natural undamped frequency by./| 2, and decreases 
the damping ratio by p 

The first condition can be seen from the following considerations. 
From Chapter II, we have shown that one version of the performance 


equation of tne system is 


l(rot) А + Са А + 5° (ез) [е;3] 5(®а) [12;м] 5^ (er) [ele] "(sg) [А;е] τ 


- 


2 
= ЖЕЛІ S (en) RE E a 6-28 


Thus, in the steady state, 


а 
3 c 6-29 
(sr) [e250] (sg) [Ase] 


A(max) = 


Since the square root sensitivity will consist of a linear 
а-с amplifier demodulator combination in cascade with a passive attenu- 
ating network which will, in the ideal case, follow a square root functicn. 


Hence, we may write 
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S tex [else] | um aie (att) [eb se] 6-30 


апа 5 
RSS 
ibm. Ln о lugo] ο -υ 6-31 


A(max) 
max 5^ (amp) [ee] S" (att) [ese] °(sg) [Ase | 


But, since 


S(ovss)|A;u | *  (tgz)|i2;w] р “η 9 Gto eise] "coo [je] 


6-32 
A(max) 7 E^ (max) Stg) lim] 5" (ea) le; 5-33 
S(ovss) | 4;M | 
or finally, 
(max) = 5 const, £m 
(ovss) [45M] 


Thus, we have the two results that increasing the sensitivity 
of the amplifier-demodulator in cascade with the non-linear attenuator 
by the factor X 2, doubles the loop gain and cuts the shaft angle in 
half. 

Referring again to Fig. 6-3, we have shown thet at any poirt 


in the (A, Eg) plane 


ОБО во 6-35 
дА y f(A) 
and 
dE, - 22,600 220 
ЧА | f(A) 7 


in order to stay within our 25% tolerances on the natural frequency and 


damping ratio of the system. 





Starting at the point, (109, 70.7V), 6-85is represented by 
curve | and 6-6 139 represented by curve 3. It will be easier to design 
a non-linear device if we plot 6-Sand 6<%as shown in Fig. 6-4. Аз 
long as we stay within the curves (A) and (B) as E, is changed, then 
we are meeting the required tolerances on W, and (DR) (ovss) - In the 
laboratory test system, we shall use a straight line attenuator which 
will abruptly bend at points where the slope of the original line will 
cause the tolerances on the dynamics to be exceeded if carried any 
further. 

Thus, starting at the point, (70.7V, 320V/rad), on Fig. 6-4, 
we proceed along the constant slope III until E, is lowered to 25.5V 
where we must break in order to stay within the required boundary 
curves À and B. The same procedure is carried out for slopes II and 
I. Slope I is carried straight through to the Eo = O line and the 
dynamic tolerances are exceeded at Eo = 3.2V. Because of the infinite 
slope required at the origin, the best we can do is to make the slope 
for low signal levels as steep as possible and specify for what levels 
the tolerances on Wy and (DR) (ovss) are exceeded. Fig. 6-3 shows how 
curves I, II, and III look on the A, f(A) plane. An important point 
to be brought out here is that slope I crosses the abscissa at some 
small positive rotation of A. Thus, we have achieved our goal of 
putting a stop (electrical in this case) so as to keep A from passing 
into the unstable region. Our next step is to instrumentize a suitable 
system to give the response represented by slopes 1, 11, and III. 

Our first step in developing an actual circuit that will 
represent the three different regions is to determine from Fig. 6-3 


the actual gains that are required for each region. What will be done 


89 


FIGURE 6-4 
VS. Eo FOR THE GIVEN 


dE 
dA 





<> 
P d 
«Ὁ 
D 
шы 
= 
ча 
ows 
L.J 
= 
== 
< 
v» 
eS 
= 
ос 
шш 
—Jj 
со 
j= 
< 
S 
ceo 


| у 


PLOT OF 


TRENT 
ШІ! 1 
ΗΕ. 
ΠΠ... ΠΗ 
ү Шш] 
ШІСІ; 

ΠΗ 


Ο 
———- 


Eo = f (A) 


CL 

ЦВЕ 
ерер T | 
Ре 
о Сис ЕЕЕ ЕЗИ a T aa Bc TN Ear етае ee иы АЩ, 5) 

τ ΓΕ! ΕΕ. 


= 


сы 
РИ 


1% 


= = 
— 
Б 
a 
ТА 
ы 
ІНІ 
В = 
rers 
ας. 
m 
Е 
J 
ж 
Е 


НТ ШШЕ 


-4— (V) }j = > 


10,000 74 


will be to attenuate the amplifier guin required for step I when 
passing into regions II and III. 
Taking the slopes from Fig. 6-3 which will give us the full 


ten degree of shaft rotation for maximum Eg we have 


Slope 5 (атр) Е;е 


T 192.3 v/v 
[Т 68.1 ν/ν 
ТЇЇ 23.9 у/у 


If we desire five degrees for maximum shaft angle, then each 
of the above gains are multiplied by the square root of two. 

The experimental results to follow will also show how the 
dynamics vary quite radically when a purely linear device is used in 
the feedback loop. For ten degree operation this characteristic is 
represented by slope 6 of Fig. 6-3 and C of Fig. 6-4. Frim Fig. 6-4, 
we see that the required 25% tolerances are met for values of Eo 
ranging from 20V to 56V. For values of Eo greater than 56V, the 
slope of the linear feedback device is too high as can be seen from 
Fig. 6-4 and hence the natural undamped frequency is more than 25% 
greater than Wo (ref) and conversely the damping ratio becomes less than 
ΠΕ сузе) (тег) - 25% (DR) ovss) (ref) ‘ Just the opposite is true for 
О ЕЕ $ 20V. 

Fig. 6-5 shows the attenuator used in the laboratory test 
system which, in conjunction with a conventional a-c amplifier and 
demodulator, will follow slopes I, II, and III of Fig. 6-3. 

In order to design the attenuator of Fig. 6-5, we need the 
equations of slopes I, II, and III, which are obtained as follows from 


Fig. 6-3. 
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LABORATORY NON-LINEAR ATTENUATOR 
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Slope I : Eg (апр) - 29 


Slope II : KE, 0354 E amp) - 9.91 


Slope III: E, = „125 E( amp) + 14.7 


Since the maximum slope to be obtained from the attenuator shown in 
Fig. 6-5 is .5, we shall double the original value for the gain of 
slope I to compensate for this. Thus, the actual equations to be 
fitted by the attenuator are 
Slope I : Е 
Slope II : Eo 


Slope III: Ep 


The diode switching circuit shown in Fig. 6-5 will follow 
the double line curve shown in Fig. 6-7. This characteristic plus the 
negative 14,.5 volts is then summed by the feedback amplifier to give 
the desired characteristic. In the actual circuit, all polarities are 
reversed since the DC amplifier involves a 180 degree phase change and 
only positive outputs of the DC amplifier are possible because of the 
unidirectional flow of current in the current source vacuum tube. 

Specifically, the operation of the diode switching circuit is 
explained as follows. As the input, E(amp) » increases from zero along 
slope I, neither diode is conducting and the voltage at A fed into the 


summing amplifier is 


Ед = Ep E(amp) ^ *? B(amp) 


However, when A= ,0175 rad as seen from Fig. 6-7, it is 


necessary to change this characteristic to follow slope II. Thus, at 


93 


| 70.7 
5 | 
= = 
iue E ға SLOPE = .124 
SLOPE = .354 | 
8.5 
4 | 
z Z SLOPE = 5 A (rot) (rad) ‚174 
/ 76.1 58.8 132 4 
| =| (amp) (volts) 
-29 | 


FIGURE 6-6 


DESIRED CHARACTERISTIC TO BE 
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dlp 


this point, V-l goes into operation and the slope is now determined to 


a very close approximation by 


(Ro + На) (Б) 
o EHE M = BT 


К] + (Ro + Ез) (6) 
Ro + R3 + R, 


In the above expression, we have neglected the small resistance of the 
6AL5 diode (varying from 200 to 1000 ohms) and that of the biasing 
network. 
Going through elementary calculations it is found that when 
А = .0175 гаа, the desired switching point, then the voltage at point A 
is - 7.3V. Thus, with the plate of V-1 set at ~ 7.3V as shown, the 
correct breakpoint is determined. Slope II is followed until A = .0335 rad. 


for the full ten degree operation. At this point, we desire V-2 to 





start conducting thus effectively shorting δη. Тре slope for ËA 
Е(ашр) 
is then quite closely found to be 


Ro Ry, 
m oS; 200 
πι + Ro Ry, 
Ro + Ry 

Again, elementary calculations have shown that when A= .0335, 
point B of Fig. 6-5 is at - 9V. Thus, setting the plate at - 9V will 
give the proper switching action. 

Fig. 6-8 shows the response of the orthogonal vector summing 
system to 3 volt step functions without a non-linear attenuator in the 
feedback loop. А Sanborn pen recorder was used to plot the output of 


the system as a function of time in seconds. The gain was set so that 
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70.7 volts output gave the full ten degree operation. Аз explained 
previously, plot C of Fig. 6-4 shows the straight gain characteristic 
for full ten degree operation. Only at E, = 34 volts is the slope of 
the linear characteristic the same as that of the ideal curve. This 
is the point on Fig. 6-4 where dj - do since the ideal curve is half 
way between A and B. 

From Fig. 6-8, it appears that the r oss referring to 
the linear second order system at E, - 34V is about .6. 


To check this, we make use of the theoretical expression 


5 а БЕ ҒА 
(DR) (ovsa) 2 5(оузз) [А;м] 1 τ 


Cg has been found to be 2.5 x 10% с at 2009F. This 
rad/sec 


was calculated by putting a known small current through a torque 
generator winding and measuring the time taken το travel through a 
given angle. 

Substituting numerical values into 6-37 gives 


(DR) .625 which is a good check on the experimental 


ovss) ` 
result. 
Wn, has already been found to be 60 rad/sec for the ideal case. 


Thus, a rough check on the time it takes for the output to first pass 


through its final asymptotic value can be found by the expression. 


27 
Тр = & х = 
ог 
т 
ТЕ — 6-38 


aw, «| 1 - (DR)? 


In the case of the ideal systen, Тр = ,033 sec, a value that, 


98 


should correspond to the 34V point.  Interpolating between the 30V 
and 40V points on Fig. 6-8, this value is fairly well checked. 

Fig. 6-9 shows the same straight gain type of system as that 
of Fig. 6-8 except that Éo(max) corresponds to 5° of shaft angle instead 
of the full 10° in the previous case. Thus, at corresponding input 
levels, the damping ratio associated with the linear squared system of 


i 
N2 


damped natural frequencies of the levels of Fig. 6-9 are Ἡ 2 times 





Fig. 6-9 is times those of Fig. 6-8. Оп the other hand, the un- 
those of the corresponding levels of Fig. 6-8. A&A good checkpoint for 
copparing the two systems is at the 30V point. From Fig. 6-8, the 
damping ratio can be seen to be about „7. Оп Fig. 6-9, it appears that 


the damping ratio of the same level is between .4 and .§. But —}- 


2 
times .7 is .49. Because of the small scale required of the recorder, 


the data is necessarily rough and thus it is quite proper to speak of 
any constants derived from the experimental results only in approximate 
terms. 

We shall now use the non-linear attenuator which was described 
earlier in this chapter. Figs. 6-104 and 6-10B show the responses of 
the system to a 3V step function for various operating levels. It is 
to be noted that for all levels from 3V to 50V, the dynamics are 
essentially constant. Even at the zero volt level, the response is 
greatly improved due to the very high gain at the very low levels. | 

Maxima were limited to 50V in the experimental results since 
only one input was used and each input channel was designed to take a 
maximum of 50V. 

In conclusion, it is hoped that a fairly logical approach has 


been made in analyzing this particular non-linear system. Admittedly 
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the title, "A Study of the Dynamics Associated With A Non-linear 
Closed Loop System," left a certain amount of leeway in just what 
was to be accomplished. However, the instrumentation of the ideal 
system in Chapter II in addition to the development of a method for 
designing a practical system shown in Chapter III should cover most 
of the background necessary to develop a practical system similar to 
that shown in Chapter VI. 

The author developed a firm appreciation for the difficulties 
involved in analyzing non-linear systems, and why they are avoided if 
possible. Although the use of perturbation theory is the most general 
method for attacking a non-linear system, eacn non-linear system is a 
separate problem in itself and must be treated as such. Thus, the 
method used to attack the Orthogonal Vector Summing System using 
microsyns probably will apply to very few if any, other non-linear 
systems. 

Аз can be seen from the data sheets, the results from the 
laboratory test system seemed to agree very well with the gereral 
theory developed prior to Chapter VI. Of course, there are probably 
better and easier methods for obtaining the required non-linear attenu- 
ator than the diode switching circuits. A light bulb could very well 
have the desired non-linear characteristics so that the pre-assigned 
tolerances on the dynamics will be maintained throughout the operating 
range, (except near zero). However, the diode switching circuit does 
not take much space and is fairly easy to design. Furthermore, by 
using more diode switches, the ideal square root curve can be approached 


as closely as desired. 
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APPENDIX 


Brief Explanation of Terminology 
Used in This Thesis 


Given: 1. Ап operating component, oc. 
2. An input quantity, q;. 
3. An output quantity, qo. 


Then the operating component performance function, (PF) (oc); 

relates the input and output quantities as follows: 
do = (РР)(ос) 9 

(PF) (oc) can be separated into a static sensitivity, 
*(oc) [31590] » which is a constant of proportionality relating the input 
&nd output quantities, and a frequency function, (ЕР) (ос), which tells 
how до is effected dynamically by the operating component. 

For example, if we have a torque generator acting on a shaft, 
referred here as a torque summing member, (tsm), then 

Μεσα) = S(tg) [325m] 2°(in) - 

This equation is read, torque applied to the torque summing member is 
equal to the product of the sensitivity of the torque generator for input 
current squared and torque out multiplied by the square of the input 
current, 


The frequency function will be of the form 


а. ра 21 №... 1 
лы ЕР... 





к 
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where p, as used in this thesis, is the Laplace transform variable. 
For further explanation of the Laplace transform, the reader is referred 


to Reference 4 of the Bibliography, 
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